"MODELO WRIGHT-FISHER DE LA DERIVA GENÉTICA".

La matriz que define el modelo sería:

A:=(({{1, 0.316, 0.062, 0.004, 0}, {0, 0.422, 0.250, 0.047, 0}, {0, 0.211, 0.376, 0.211, 0}, {0, 0.047, 0.250, 0.422, 0}, {0, 0.004, 0.062, 0.316, 1}}))

MatrixForm[MatrixPower[A, 6]]

La cadena de Markov no es regular. Podemos calcular la matriz potencia diagonalizando la matriz, (A^k)=C.(D^k).(C^(-1))

Eigenvalues[A]

Eigenvectors[A]

paso := Transpose[Eigenvectors[A]]

MatrixForm[paso]

diagonal:=(({{1, 0, 0, 0, 0}, {0, 1, 0, 0, 0}, {0, 0, 0.75019, 0, 0}, {0, 0, 0, 0.375, 0}, {0, 0, 0, 0, 0.0938}}))

potencia := 
      paso . DiagonalMatrix[{1, 1, 0.75019^100, 0.375^100, 0.0938^100}] . 
        Inverse[paso]

MatrixForm[potencia]

O bien diréctamente,

MatrixForm[MatrixPower[A, 100]]

Como inicialmente el número de alelos es de dos, el vector inicial será:

(Subscript[X, 0]):=(({{0}, {0}, {1}, {0}, {0}}))

Después de una generación:

MatrixForm[A . (Subscript[X, 0])]

Después de 100 generaciones

MatrixForm[MatrixPower[A, 100] . (Subscript[X, 0])]

Resultado: 
La tendencia de la población cuando pasen un número grande de generaciones será  del 50 % para ningún alelos A1 y otro 50 % para cuatro alelos A1. Es decir, fijar la homocigosis y perder la heterocigosis.